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FORMATION OF A FIELD REVERSED CONFIGURATION FOR MAGNETIC AND 
ELECTROSTATIC CONFINEMENT OF PLASMA 

This invention was made with Government support under Contract No. N00014-99-1- 

0857, awarded by the Office of Naval Research. Some background research was supported 
5 by the U.S. Department of Energy for 1992 to 1993. The Government has certain rights in 

this invention. 

CROSS-REFERENCE TO RELATED APPLICATIONS 

This application is a continuation of U.S. Serial No. 10/328,701 filed December 23, 
2002, which is a continuation of U.S. Serial No. 10/066,424, filed January 31, 2002, which 
10 applications are fully incorporated herein by reference, and which is related to provisional 
U.S. application serial no. 60/266,074, filed February 1, 2001, and provisional U.S. 
application serial no. 60/297,086, filed on June 8, 2001, which applications are fully 
incorporated herein by reference. 

FIELD OF THE INVENTION 

1 5 The invention relates generally to the field of plasma physics, and, in particular, to 

methods and apparati for confining plasma. Plasma confinement is particularly of interest for 
the purpose of enabling a nuclear fusion reaction. 
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BACKGROUND OF THE INVENTION 

Fusion is the process by which two light nuclei combine to form a heavier one. The 
fusion process releases a tremendous amount of energy in the form of fast moving particles. 
Because atomic nuclei are positively charged — due to the protons contained therein — there 
is a repulsive electrostatic, or Coulomb, force between them. For two nuclei to fuse, this 
repulsive barrier must be overcome, which occurs when two nuclei are brought close enough 
together where the short-range nuclear forces become strong enough to overcome the 
Coulomb force and fuse the nuclei. The energy necessary for the nuclei to overcome the 
Coulomb barrier is provided by their thermal energies, which must be very high. For 
example, the fusion rate can be appreciable if the temperature is at least of the order of 10 4 
eV — corresponding roughly to 100 million degrees Kelvin. The rate of a fusion reaction is a 
function of the temperature, and it is characterized by a quantity called reactivity. The 
reactivity of a D-T reaction, for example, has a broad peak between 30 keV and 100 keV. 
Typical fusion reactions include: 

D + D -> He 3 (0.8 MeV) + «(2.5 MeV), 
D + T a(3.6 MeV) + «(14.1 MeV), 
D + He 3 -» a(3.7 MeV) + p(UJ MeV), and 
p + B n 3<x(8.7 MeV), 
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where D indicates deuterium, T indicates tritium, a indicates a helium nucleus, n indicates a 
neutron, p indicates a proton, He indicates helium, and B 1 1 indicates Boron- 1 1 . The numbers 
in parentheses in each equation indicate the kinetic energy of the fusion products. 

The first two reactions listed above — the D-D and D-T reactions — are neutronic, 
5 which means that most of the energy of their fusion products is carried by fast neutrons. The 

disadvantages of neutronic reactions are that (1) the flux of fast neutrons creates many 
problems, including structural damage of the reactor walls and high levels of radioactivity for 
most construction materials; and (2) the energy of fast neutrons is collected by converting 
their thermal energy to electric energy, which is very inefficient (less than 30%). The 

10 advantages of neutronic reactions are that (1) their reactivity peaks at a relatively low 

temperature; and (2) their losses due to radiation are relatively low because the atomic 
numbers of deuterium and tritium are 1. 

The reactants in the other two equations — D-He 3 and p-B n — are called advanced 
fuels. Instead of producing fast neutrons, as in the neutronic reactions, their fusion products 

15 are charged particles. One advantage of the advanced fuels is that they create much fewer 

neutrons and therefore suffer less from the disadvantages associated with them. In the case of 
D-He 3 , some fast neutrons are produced by secondary reactions, but these neutrons account 
for only about 10 per cent of the energy of the fusion products. The p-B u reaction is free of 
fast neutrons, although it does produce some slow neutrons that result from secondary 

20 reactions but create much fewer problems. Another advantage of the advanced fuels is that 
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the energy of their fusion products can be collected with a high efficiency, up to 90 per cent. 
In a direct energy conversion process, their charged fusion products can be slowed down and 
their kinetic energy converted directly to electricity. 

The advanced fuels have disadvantages, too. For example, the atomic numbers of the 
advanced fuels are higher (2 for He 3 and 5 for B 11 ). Therefore, their radiation losses are 
greater than in the neutronic reactions. Also, it is much more difficult to cause the advanced 
fuels to fuse. Their peak reactivities occur at much higher temperatures and do not reach as 
high as the reactivity for D-T. Causing a fusion reaction with the advanced fuels thus 
requires that they be brought to a higher energy state where their reactivity is significant. 
Accordingly, the advanced fuels must be contained for a longer time period wherein they can 
be brought to appropriate fusion conditions. 

The containment time for a plasma is At = r 2 /D, where r is a minimum plasma 
dimension and D is a diffusion coefficient. The classical value of the diffusion coefficient is 
D c = a] '/ x /e , where a t is the ion gyroradius and x ie is the ion-electron collision time. 
Diffusion according to the classical diffusion coefficient is called classical transport. The 
Bohm diffusion coefficient, attributed to short-wavelength instabilities, is D B = (l/16)a?Q., 
where Q, is the ion gyro frequency. Diffusion according to this relationship is called 
anomalous transport. For fusion conditions, D B I D c = (l/16)Q,x /e = 10 8 , anomalous transport 
results in a much shorter containment time than does classical transport. This relation 
determines how large a plasma must be in a fusion reactor, by the requirement that the 
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containment time for a given amount of plasma must be longer than the time for the plasma 
to have a nuclear fusion reaction. Therefore, classical transport condition is more desirable in 
a fusion reactor, allowing for smaller initial plasmas. 

In early experiments with toroidal confinement of plasma, a containment time of At = 
was observed. Progress in the last 40 years has increased the containment time to At = 
1000 r^lD B . One existing fusion reactor concept is the Tokamak. The magnetic field of a 
Tokamak 68 and a typical particle orbit 66 are illustrated in Fig. 5. For the past 30 years, 
fusion efforts have been focussed on the Tokamak reactor using a D-T fuel. These efforts 
have culminated in the International Thermonuclear Experimental Reactor (ITER), illustrated 

10 in Fig. 7. Recent experiments with Tokamaks suggest that classical transport, At = r 2 /D c , is 

possible, in which case the minimum plasma dimension can be reduced from meters to 
centimeters. These experiments involved the injection of energetic beams (50 to 100 keV), to 
heat the plasma to temperatures of 10 to 30 keV. See W. Heidbrink & G. J. Sadler, 34 
Nuclear Fusion 535 (1994). The energetic beam ions in these experiments were observed to 

15 slow down and diffuse classically while the thermal plasma continued to diffuse anomalously 

fast. The reason for this is that the energetic beam ions have a large gyroradius and, as such, 
are insensitive to fluctuations with wavelengths shorter than the ion gyroradius (k < ai). The 
short-wavelength fluctuations tend to average over a cycle and thus cancel. Electrons, 
however, have a much smaller gyroradius, so they respond to the fluctuations and transport 

20 anomalously. 
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Because of anomalous transport, the minimum dimension of the plasma must be at 
least 2.8 meters. Due to this dimension, the ITER was created 30 meters high and 30 meters 
in diameter. This is the smallest D-T Tokamak-type reactor that is feasible. For advanced 
fuels, such as D-He 3 and p-B { \ the Tokamak-type reactor would have to be much larger 
because the time for a fuel ion to have a nuclear reaction is much longer. A Tokamak reactor 
using D-T fuel has the additional problem that most of the energy of the fusion products 
energy is carried by 14 MeV neutrons, which cause radiation damage and induce reactivity in 
almost all construction materials due to the neutron flux. In addition, the conversion of their 
energy into electricity must be by a thermal process, which is not more than 30% efficient. 

Another proposed reactor configuration is a colliding beam reactor. In a colliding 
beam reactor, a background plasma is bombarded by beams of ions. The beams comprise 
ions with an energy that is much larger than the thermal plasma. Producing useful fusion 
reactions in this type of reactor has been infeasible because the background plasma slows 
down the ion beams. Various proposals have been made to reduce this problem and 
maximize the number of nuclear reactions. 

For example, U.S. Patent No. 4,065,351 to Jassby et al. discloses a method of 
producing counterstreaming colliding beams of deuterons and tritons in a toroidal 
confinement system. In U.S. Patent No. 4,057,462 to Jassby et al., electromagnetic energy is 
injected to counteract the effects of bulk equilibrium plasma drag on one of the ion species. 
The toroidal confinement system is identified as a Tokamak. In U.S. Patent No. 4,894,199 to 
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Rostoker, beams of deuterium and tritium are injected and trapped with the same average 
velocity in a Tokamak, mirror, or field reversed configuration. There is a low density cool 
background plasma for the sole purpose of trapping the beams. The beams react because they 
have a high temperature, and slowing down is mainly caused by electrons that accompany the 
5 injected ions. The electrons are heated by the ions in which case the slowing down is 

minimal. 

In none of these devices, however, does an equilibrium electric field play any part. 
Further, there is no attempt to reduce, or even consider, anomalous transport. 

Other patents consider electrostatic confinement of ions and, in some cases, magnetic 

10 confinement of electrons. These include U.S. Patent No. 3,258,402 to Farnsworth and U.S. 
Patent No. 3,386,883 to Farnsworth, which disclose electrostatic confinement of ions and 
inertial confinement of electrons; U.S. Patent No. 3,530,036 to Hirsch et al. and U.S. Patent 
No. 3,530,497 to Hirsch et al. are similar to Farnsworth; U.S. Patent No. 4,233,537 to 
Limpaecher, which discloses electrostatic confinement of ions and magnetic confinement of 

15 electrons with multipole cusp reflecting walls; and U.S. Patent No. 4,826,646 to Bussard, 

which is similar to Limpaecher and involves point cusps. None of these patents consider 
electrostatic confinement of electrons and magnetic confinement of ions. Although there 
have been many research projects on electrostatic confinement of ions, none of them have 
succeeded in establishing the required electrostatic fields when the ions have the required 
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density for a fusion reactor. Lastly, none of the patents cited above discuss a field reversed 
configuration magnetic topology. 

The field reversed configuration (FRC) was discovered accidentally around 1960 at 
the Naval Research Laboratory during theta pinch experiments. A typical FRC topology, 
wherein the internal magnetic field reverses direction, is illustrated in Fig. 8 and Fig. 10, and 
particle orbits in a FRC are shown in Fig. 1 1 and Fig. 14. Regarding the FRC, many research 
programs have been supported in the United States and Japan. There is a comprehensive 
review paper on the theory and experiments of FRC research from 1960-1988. See M. 
Tuszewski, 28 Nuclear Fusion 2033, (1988). A white paper on FRC development describes 
the research in 1996 and recommendations for future research. See L. C. Steinhauer et al., 30 
Fusion Technology 116 (1996). To this date, in FRC experiments the FRC has been formed 
with the theta pinch method. A consequence of this formation method is that the ions and 
electrons each carry half the current, which results in a negligible electrostatic field in the 
plasma and no electrostatic confinement. The ions and electrons in these FRCs were 
contained magnetically. In almost all FRC experiments, anomalous transport has been 
assumed. See, e.g., Tuszewski, beginning of section 1.5.2, at page 2072. 

SUMMARY OF THE INVENTION 

To address the problems faced by previous plasma containment systems, a system and 
apparatus for containing plasma are herein described in which plasma ions are contained 
magnetically in stable, large orbits and electrons are contained electrostatically in an energy 
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well. A major innovation of the present invention over all previous work with FRCs is the 
simultaneous electrostatic confinement of electrons and magnetic confinement of ions, which 
tends to avoid anomalous transport and facilitate classical containment of both electrons and 
ions. In this configuration, ions may have adequate density and temperature so that upon 
collisions they are fused together by the nuclear force, thus releasing fusion energy. 

In a preferred embodiment, a plasma confinement system comprises a chamber, a 
magnetic field generator for applying a magnetic field in a direction substantially along a 
principle axis, and an annular plasma layer that comprises a circulating beam of ions. Ions of 
the annular plasma beam layer are substantially contained within the chamber magnetically in 
orbits and the electrons are substantially contained in an electrostatic energy well. In one 
aspect of one preferred embodiment a magnetic field generator comprises a current coil. 
Preferably, the system further comprises mirror coils near the ends of the chamber that 
increase the magnitude of the applied magnetic field at the ends of the chamber. The system 
may also comprise a beam injector for injecting a neutralized ion beam into the applied 
magnetic field, wherein the beam enters an orbit due to the force caused by the applied 
magnetic field. In another aspect of the preferred embodiments, the system forms a magnetic 
field having a topology of a field reversed configuration. 

Also disclosed is a method of confining plasma comprising the steps of magnetically 
confining the ions in orbits within a magnetic field and electrostatically confining the 
electrons in an energy well. An applied magnetic field may be tuned to produce and control 
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the electrostatic field. In one aspect of the method the field is tuned so that the average 
electron velocity is approximately zero. In another aspect, the field is timed so that the 
average electron velocity is in the same direction as the average ion velocity. In another 
aspect of the method, the method forms a field reversed configuration magnetic field, in 
which the plasma is confined. 

In another aspect of the preferred embodiments, an annular plasma layer is contained 
within a field reversed configuration magnetic field. The plasma layer comprises positively 
charged ions, wherein substantially all of the ions are non-adiabatic, and electrons contained 
within an electrostatic energy well. The plasma layer is caused to rotate and form a magnetic 
self-field of sufficient magnitude to cause field reversal. 

In other aspects of the preferred embodiments, the plasma may comprise at least two 
different ion species, one or both of which may comprise advanced fuels. 

Having a non-adiabatic plasma of energetic, large-orbit ions tends to prevent the 
anomalous transport of ions. This can be done in a FRC, because the magnetic field vanishes 
(i.e., is zero) over a surface within the plasma. Ions having a large orbit tend to be insensitive 
to short-wavelength fluctuations that cause anomalous transport. 

Magnetic confinement is ineffective for electrons because they have a small 
gyroradius — due to their small mass — and are therefore sensitive to short-wavelength 
fluctuations that cause anomalous transport. Therefore, the electrons are effectively confined 
in a deep potential well by an electrostatic field, which tends to prevent the anomalous 
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transport of energy by electrons. The electrons that escape confinement must travel from the 
high density region near the null surface to the surface of the plasma. In so doing, most of 
their energy is spent in ascending the energy well. When electrons reach the plasma surface 
and leave with fusion product ions, they have little energy left to transport. The strong 
electrostatic field also tends to make all the ion drift orbits rotate in the diamagnetic direction, 
so that they are contained. The electrostatic field further provides a cooling mechanism for 
electrons, which reduces their radiation losses. 

The increased containment ability allows for the use of advanced fuels such as D-He 3 
and p-B u 9 as well as neutronic reactants such as D-D and D-T. In the D-He 3 reaction, fast 
neutrons are produced by secondary reactions, but are an improvement over the D-T reaction. 
The p-B u reaction, and the like, is preferable because it avoids the problems of fast neutrons 
completely. 

Another advantage of the advanced fuels is the direct energy conversion of energy 
from the fusion reaction because the fusion products are moving charged particles, which 
create an electrical current. This is a significant improvement over Tokamaks, for example, 
where a thermal conversion process is used to convert the kinetic energy of fast neutrons into 
electricity. The efficiency of a thermal conversion process is lower than 30%, whereas the 
efficiency of direct energy conversion can be as high as 90%. 
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Other aspects and features of the present invention will become apparent from 
consideration of the following description taken in conjunction with the accompanying 
drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Preferred embodiments are illustrated by way of example, and not by way of 
limitation, in the figures of the accompanying drawings, in which like reference numerals 
refer to like components. 

Figs. 1 A and IB show, respectively, the Lorentz force acting on a positive and a 
negative charge. 

Figs. 2A and 2B show Larmor orbits of charged particles in a constant magnetic field. 

Fig. 3 shows the ExB drift. 

Fig. 4 shows the gradient drift. 

Fig. 5 shows an adiabatic particle orbit in a Tokamak. 

Fig. 6 shows a non-adiabatic particle orbit in a betatron. 

Fig. 7 shows the International Thermonuclear Experimental Reactor (ITER). 

Fig. 8 shows the magnetic field of a FRC. 

Figs. 9A and 9B show, respectively, the diamagnetic and the counterdiamagnetic 
direction in a FRC. 

Fig. 10 shows the colliding beam system. 
Fig. 1 1 shows a betatron orbit. 
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Figs. 12 A and 12B show, respectively, the magnetic field and the direction of the 
gradient drift in a FRC. 

Figs. 13A and 13B show, respectively, the electric field and the direction of the ExB 
drift in a FRC. 

Figs. 14 A, 14B and 14C show ion drift orbits. 

Figs. 15A and 15B show the Lorentz force at the ends of a FRC. 

Figs. 16A and 16B show the tuning of the electric field and the electric potential in 
the colliding beam system. 

Fig. 17 shows a Maxwell distribution. 

Figs. 18A and 18B show transitions from betatron orbits to drift orbits due to large- 
angle, ion-ion collisions. 

Figs. 19 show A, B, C and D betatron orbits when small-angle, electron-ion collisions 
are considered. 

Figs. 20A, 20B and 20C show the reversal of the magnetic field in a FRC. 
Figs. 21 A, 21B, 21C and 21D show the effects due to tuning of the external magnetic 
field Bo in a FRC. 

Figs. 22A, 22B, 22C and 22D show iteration results for a D-T plasma. 
Figs. 23A, 23B, 23C, and 23D show iteration results for a D-He 3 plasma. 
Fig. 24 shows iteration results for a p-B u plasma. 
Fig. 25 shows an exemplary confinement chamber. 
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Fig. 26 shows a neutralized ion beam as it is electrically polarized before entering a 
confining chamber. 

Fig. 27 is a head-on view of a neutralized ion beam as it contacts plasma in a 
confining chamber. 

Fig. 28 is a side view schematic of a confining chamber according to a preferred 
embodiment of a start-up procedure. 

Fig. 29 is a side view schematic of a confining chamber according to another 
preferred embodiment of a start-up procedure. 

Fig. 30 shows traces of B-dot probe indicating the formation of a FRC. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

An ideal fusion reactor solves the problem of anomalous transport for both ions and 
electrons. The anomalous transport of ions is avoided by magnetic confinement in a field 
reversed configuration (FRC) in such a way that the majority of the ions have large, non- 
adiabatic orbits, making them insensitive to short-wavelength fluctuations that cause 
anomalous transport of adiabatic ions. For electrons, the anomalous transport of energy is 
avoided by tuning the externally applied magnetic field to develop a strong electric field, 
which confines them electrostatically in a deep potential well. Moreover, the fusion fuel 
plasmas that can be used with the present confinement process and apparatus are not limited 
to neutronic fuels only, but also advantageously include advanced fuels. (For a discussion of 
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advanced fuels, see R. Feldbacher & M. Heindler, Nuclear Instruments and Methods in 
Physics Research, A271(1988)JJ-64 (North Holland Amsterdam).) 

The solution to the problem of anomalous transport found herein makes use of a 
specific magnetic field configuration, which is the FRC. In particular, the existence of a 
region in a FRC where the magnetic field vanishes makes it possible to have a plasma 
comprising a majority of non-adiabatic ions. 
Background Theory 

Before describing the system and apparatus in detail, it will be helpful to first review 
a few key concepts necessary to understand the concepts contained herein. 
Lorentz Force and Particle Orbits in a Magnetic Field 

A particle with electric charge q moving with velocity v in a magnetic field B 
experiences a force F L given by 

c 

The force F L is called the Lorentz force. It, as well as all the formulas used in the present 
discussion, is given in the gaussian system of units. The direction of the Lorentz force 
depends on the sign of the electric charge q. The force is perpendicular to both velocity and 
magnetic field. Fig. 1 A shows the Lorentz force 30 acting on a positive charge. The velocity 
of the particle is shown by the vector 32. The magnetic field is 34. Similarly, Fig. IB shows 
the Lorentz force 30 acting on a negative charge. 
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As explained, the Lorentz force is perpendicular to the velocity of a particle; thus, a 
magnetic field is unable to exert force in the direction of the particle's velocity. It follows 

from Newton's second law, F = ma y that a magnetic field is unable to accelerate a particle in 
the direction of its velocity. A magnetic field can only bend the orbit of a particle, but the 
magnitude of its velocity is not affected by a magnetic field. 

Fig. 2A shows the orbit of a positively charged particle in a constant magnetic field 
34. The Lorentz force 30 in this case is constant in magnitude, and the orbit 36 of the particle 
forms a circle. This circular orbit 36 is called a Larmor orbit. The radius of the circular orbit 
36 is called a gyroradius 38. 

Usually, the velocity of a particle has a component that is parallel to the magnetic 
field and a component that is perpendicular to the field. In such a case, the particle undergoes 
two simultaneous motions: a rotation around the magnetic field line and a translation along 
it. The combination of these two motions creates a helix that follows the magnetic field line 
40. This is indicated in Fig. 2B. 

A particle in its Larmor orbit revolves around a magnetic field line. The number of 
radians traveled per unit time is the particle's gyro frequency, which is denoted by Q and 
given by 

a-**. (2) 

mc 
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where m is the mass of the particle and c is the speed of light. The gyroradius a L of a charged 
particle is given by 

(3) 

where vi is the component of the velocity of the particle perpendicular to the magnetic field. 
ExB Drift and Gradient Drift 

Electric fields affect the orbits of charged particles, as shown in Fig. 3. In Fig. 3, the 
magnetic field 44 points toward the reader. The orbit of a positively charged ion due to the 
magnetic field 44 alone would be a circle 36; the same is true for an electron 42. In the 
presence of an electric field 46, however, when the ion moves in the direction of the electric 
field 46, its velocity increases. As can be appreciated, the ion is accelerated by the force qE . 
It can further be seen that, according to Eq. 3, the ion's gyroradius will increase as its 
velocity does. 

As the ion is accelerated by the electric field 46, the magnetic field 44 bends the ion's 
orbit. At a certain point the ion reverses direction and begins to move in a direction opposite 
to the electric field 46. When this happens, the ion is decelerated, and its gyroradius 
therefore decreases. The ion's gyroradius thus increases and decreases in alternation, which 
gives rise to a sideways drift of the ion orbit 48 in the direction 50 as shown in Fig. 3. This 
motion is called ExB drift. Similarly, electron orbits 52 drift in the same direction 50. 
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A similar drift can be caused by a gradient of the magnetic field 44 as illustrated in 
Fig. 4. In Fig. 4, the magnetic field 44 points towards the reader. The gradient of the 
magnetic field is in the direction 56. The increase of the magnetic field's strength is depicted 
by the denser amount of dots in the figure. 
5 From Eqs. 2 and 3, it follows that the gyroradius is inversely proportional to the 

strength of the magnetic field. When an ion moves in the direction of increasing magnetic 
field its gyroradius will decrease, because the Lorentz force increases, and vice versa. The 
ion's gyroradius thus decreases and increases in alternation, which gives rise to a sideways 
drift of the ion orbit 58 in the direction 60. This motion is called gradient drift. Electron 
10 orbits 62 drift in the opposite direction 64. 

Adiabatic and Non-Adiabatic Particles 

Most plasma comprises adiabatic particles. An adiabatic particle tightly follows the 
magnetic field lines and has a small gyroradius. Fig. 5 shows a particle orbit 66 of an 
adiabatic particle that follows tightly a magnetic field line 68. The magnetic field lines 68 
1 5 depicted are those of a Tokamak. 

A non-adiabatic particle has a large gyroradius. It does not follow the magnetic field 
lines and is usually energetic. There exist other plasmas that comprise non-adiabatic 
particles. Fig. 6 illustrates a non-adiabatic plasma for the case of a betatron. The pole pieces 
70 generate a magnetic field 72. As Fig. 6 illustrates, the particle orbits 74 do not follow the 
20 magnetic field lines 72. 
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Radiation in Plasmas 

A moving charged particle radiates electromagnetic waves. The power radiated by 
the particle is proportional to the square of the charge. The charge of an ion is Ze, where e is 
the electron charge and Z is the atomic number. Therefore, for each ion there will be Z free 
electrons that will radiate. The total power radiated by these Z electrons is proportional to the 
cube of the atomic number (Z 3 ). 

Charged Particles in a FRC 

Fig. 8 shows the magnetic field of a FRC. The system has cylindrical symmetry with 
respect to its axis 78. In the FRC, there are two regions of magnetic field lines: open 80 and 
closed 82. The surface dividing the two regions is called the separatrix 84. The FRC forms a 
cylindrical null surface 86 in which the magnetic field vanishes. In the central part 88 of the 
FRC the magnetic field does not change appreciably in the axial direction. At the ends 90, 
the magnetic field does change appreciably in the axial direction. The magnetic field along 
the center axis 78 reverses direction in the FRC, which gives rise to the term "Reversed" in 
Field Reversed Configuration (FRC). 

In Fig. 9A, the magnetic field outside of the null surface 94 is in the direction 96. The 
magnetic field inside the null surface is in the direction 98. If an ion moves in the direction 
100, the Lorentz force 30 acting on it points towards the null surface 94. This is easily 
appreciated by applying the right-hand rule. For particles moving in the direction 102, called 
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diamagnetic, the Lorentz force always points toward the null surface 94. This phenomenon 
gives rise to a particle orbit called betatron orbit, to be described below. 

Fig. 9B shows an ion moving in the direction 104, called counterdiamagnetic. The 
Lorentz force in this case points away from the null surface 94. This phenomenon gives rise 
5 to a type of orbit called a drift orbit, to be described below. The diamagnetic direction for 

ions is counterdiamagnetic for electrons, and vice versa. 

Fig. 10 shows a ring or annular layer of plasma 106 rotating in the ions' diamagnetic 
direction 102. The ring 106 is located around the null surface 86. The magnetic field 108 
created by the annular plasma layer 106, in combination with an externally applied magnetic 
10 field 1 10, forms a magnetic field having the topology of a FRC (The topology is shown in Fig 

8). 

The ion beam that forms the plasma layer 106 has a temperature; therefore, the 
velocities of the ions form a Maxwell distribution in a frame rotating at the average angular 
velocity of the ion beam. Collisions between ions of different velocities lead to fusion 
15 reactions. For this reason, the plasma beam layer 106 is called a colliding beam system. 

Fig. 1 1 shows the main type of ion orbits in a colliding beam system, called a betatron 
orbit 1 12. A betatron orbit 1 12 can be expressed as a sine wave centered on the null circle 
114. As explained above, the magnetic field on the null circle 1 14 vanishes. The plane of 
the orbit 1 12 is perpendicular to the axis 78 of the FRC. Ions in this orbit 1 12 move in their 
20 diamagnetic direction 102 from a starting point 1 16. An ion in a betatron orbit has two 
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motions: an oscillation in the radial direction (perpendicular to the null circle 114), and a 
translation along the null circle 1 14. 

Fig. 12A is a graph of the magnetic field 1 18 in a FRC. The field 1 18 is derived using 
a one-dimensional equilibrium model, to be discussed below in conjunction with the theory 
of the invention. The horizontal axis of the graph represents the distance in centimeters from 
the FRC axis 78. The magnetic field is in kilogauss. As the graph depicts, the magnetic field 
118 vanishes at the null circle radius 120. 

As shown in Fig. 12B, a particle moving near the null circle will see a gradient 126 of 
the magnetic field pointing away from the null surface 86. The magnetic field outside the 
null circle is 122, while the magnetic field inside the null circle is 124. The direction of the 

gradient drift is given by the cross product B x V5 , where VB is the gradient of the magnetic 
field; thus, it can be appreciated by applying the right-hand rule that the direction of the 
gradient drift is in the counterdiamagnetic direction, whether the ion is outside or inside the 
null circle 128. 

Fig. 13A is a graph of the electric field 130 in a FRC. The field 130 is derived using a 
one-dimensional equilibrium model, to be discussed below in conjunction with the theory of 
the invention. The horizontal axis of the graph represents the distance in centimeters from the 
FRC axis 78. The electric field is in volts/cm. As the graph depicts, the electric field 130 
vanishes close to the null circle radius 120. 
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As shown if Fig. 13B, the electric field for ions is deconfining; it points away from 
the null surface 132,134. The magnetic field, as before, is in the directions 122,124. It can 

be appreciated by applying the right-hand rule that the direction of the ExB drift is in the 
diamagnetic direction, whether the ion is outside or inside the null surface 136. 
5 Figs. 14A and 14B show another type of common orbit in a FRC, called a drift orbit 

138. Drift orbits 138 can be outside of the null surface, as shown in Fig. 14 A, or inside it, as 

shown in Fig. 14B. Drift orbits 138 rotate in the diamagnetic direction if the E x B drift 
dominates or in the counterdiamagnetic direction if the gradient drift dominates. The drift 
orbits 138 shown in Figs. 14A and 14B rotate in the diamagnetic direction 102 from starting 
10 point 116. 

A drift orbit, as shown in Fig. 14C, can be thought of as a small circle rolling over a 
relatively bigger circle. The small circle 142 spins around its axis in the sense 144. It also 
rolls over the big circle 146 in the direction 102. The point 140 will trace in space a path 
similar to 138. 

15 Figs. 15A and 15B show the direction of the Lorentz force at the ends of a FRC. In 

Fig. 15 A, an ion is shown moving in the diamagnetic direction 102 with a velocity 148 in a 
magnetic field 150. It can be appreciated by applying the right-hand rule that the Lorentz 
force 152 tends to push the ion back into the region of closed field lines. In this case, 
therefore, the Lorentz force 152 is confining for the ions. In Fig. 15B, an ion is shown 

20 moving in the counterdiamagnetic direction with a velocity 148 in a magnetic field 150. It 
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can be appreciated by applying the right-hand rule that the Lorentz force 152 tends to push 
the ion into the region of open field lines. In this case, therefore, the Lorentz force 152 is 
deconfining for the ions. 

Magnetic and Electrostatic Confinement in a FRC 

A plasma layer 106 (see Fig. 10) can be formed in a FRC by injecting energetic ion 
beams around the null surface 86 in the diamagnetic direction 102 of ions. (A detailed 
discussion of different methods of forming the FRC and plasma ring follows below.) In the 
circulating plasma layer 106, most of the ions have betatron orbits 1 12 (see Fig. 1 1), are 
energetic, and are non-adiabatic; thus, they are insensitive to short-wavelength fluctuations 
that cause anomalous transport. 

While studying a plasma layer 106 in equilibrium conditions as described above, it 
was discovered that the conservation of momentum imposes a relation between the angular 
velocity of ions 0% and the angular velocity of electrons co e . (The derivation of this relation is 
given below in conjunction with the theory of the invention.) The relation is 


co = co. 

e i 


CO 

1— 


^ ZeB 0 
, where Q 0 = - 


(4) 


In Eq. 4, Z is the ion atomic number, m, is the ion mass, e is the electron charge, Bo is the 
magnitude of the applied magnetic field, and c is the speed of light. There are three free 
parameters in this relation: the applied magnetic field Bo, the electron angular velocity co e , and 
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the ion angular velocity a*. If two of them are known, the third can be determined from Eq. 
4. 

Because the plasma layer 106 is formed by injecting ion beams into the FRC, the 
angular velocity of ions a* is determined by the injection kinetic energy of the beam Wj, 
5 which is given by 

W, =-m i V i 2 =-m i (o> i r 0 ) 2 . 

Here, Vi = ^r 0 , where V t is the injection velocity of ions, cch is the cyclotron frequency of 
ions, and r 0 is the radius of the null surface 86. The kinetic energy of electrons in the beam 
has been ignored because the electron mass m e is much smaller than the ion mass m,. 

10 For a fixed injection velocity of the beam (fixed the applied magnetic field B 0 can 

be tuned so that different values of co e are obtainable. As will be shown, tuning the external 
magnetic field B 0 also gives rise to different values of the electrostatic field inside the plasma 
layer. This feature of the invention is illustrated in Figs. 16A and 16B. Fig. 16A shows three 
plots of the electric field (in volts/cm) obtained for the same injection velocity, an = 1.35 x 

15 10 7 s" 1 , but for three different values of the applied magnetic field B 0 \ 


Plot 

Applied magnetic field (Bo) 

electron angular velocity (co e ) 

154 

Bo = 2.11 kG 

a> e = 0 

156 

£ 0 = 5.15kG 

o) e = 0.625 x 10 7 s" 1 

158 

B 0 = 15.5 kG 

(o e = 1.11 x 10 7 s' 1 
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The values of co e in the table above were determined according to Eq. 4. One can appreciate 
that o e > 0 means that Qo > a* in Eq. 4, so that electrons rotate in their counterdiamagnetic 
direction. Fig. 16B shows the electric potential (in volts) for the same set of values of B 0 and 
co e . The horizontal axis, in Figs. 16A and 16B, represents the distance from the FRC axis 78, 
5 shown in the graph in centimeters. The analytic expressions of the electric field and the 

electric potential are given below in conjunction with the theory of the invention. These 
expressions depend strongly on co e . 

The above results can be explained on simple physical grounds. When the ions rotate 
in the diamagnetic direction, the ions are confined magnetically by the Lorentz force. This 

10 was shown in Fig. 9 A. For electrons, rotating in the same direction as the ions, the Lorentz 
force is in the opposite direction, so that electrons would not be confined. The electrons 
leave the plasma and, as a result, a surplus of positive charge is created. This sets up an 
electric field that prevents other electrons from leaving the plasma. The direction and the 
magnitude of this electric field, in equilibrium, is determined by the conservation of 

15 momentum. The relevant mathematical details are given below in conjunction with the 

theory of the invention. 

The electrostatic field plays an essential role on the transport of both electrons and 
ions. Accordingly, an important aspect of this invention is that a strong electrostatic field is 
created inside the plasma layer 106, the magnitude of this electrostatic field is controlled by 

20 the value of the applied magnetic field Bo which can be easily adjusted. 
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As explained, the electrostatic field is confining for electrons if a> e > 0. As shown in 
Fig. 16B, the depth of the well can be increased by tuning the applied magnetic field 2? 0 . 
Except for a very narrow region near the null circle, the electrons always have a small 
gyroradius. Therefore, electrons respond to short-wavelength fluctuations with an 
5 anomalously fast diffusion rate. This diffusion, in fact, helps maintain the potential well once 

the fusion reaction occurs. The fusion product ions, being of much higher energy, leave the 
plasma. To maintain charge quasi-neutrality, the fusion products must pull electrons out of 
the plasma with them, mainly taking the electrons from the surface of the plasma layer. The 
density of electrons at the surface of the plasma is very low, and the electrons that leave the 
10 plasma with the fusion products must be replaced; otherwise, the potential well would 

disappear. 

Fig. 17 shows a Maxwellian distribution 162 of electrons. Only very energetic 
electrons from the tail 160 of the Maxwell distribution can reach the surface of the plasma 
and leave with fusion ions. The tail 160 of the distribution 162 is thus continuously created 

15 by electron-electron collisions in the region of high density near the null surface. The 

energetic electrons still have a small gyroradius, so that anomalous diffusion permits them to 
reach the surface fast enough to accommodate the departing fusion product ions. The 
energetic electrons lose their energy ascending the potential well and leave with very little 
energy. Although the electrons can cross the magnetic field rapidly, due to anomalous 

20 transport, anomalous energy losses tend to be avoided because little energy is transported. 
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Another consequence of the potential well is a strong cooling mechanism for 
electrons that is similar to evaporative cooling. For example, for water to evaporate, it must 
be supplied the latent heat of vaporization. This heat is supplied by the remaining liquid 
water and the surrounding medium, which then thermalize rapidly to a lower temperature 
5 faster than the heat transport processes can replace the energy. Similarly, for electrons, the 

potential well depth is equivalent to water's latent heat of vaporization. The electrons supply 
the energy required to ascend the potential well by the thermalization process that re-supplies 
the energy of the Maxwell tail so that the electrons can escape. The thermalization process 
thus results in a lower electron temperature, as it is much faster than any heating process. 

10 Because of the mass difference between electrons and protons, the energy transfer time from 

protons is about 1800 times less than the electron thermalization time. This cooling 
mechanism also reduces the radiation loss of electrons. This is particularly important for 
advanced fuels, where radiation losses are enhanced by fuel ions with atomic number Z > 1 . 
The electrostatic field also affects ion transport. The majority of particle orbits in the 

15 plasma layer 106 are betatron orbits 1 12. Large-angle collisions, that is, collisions with 

scattering angles between 90° and 180°, can change a betatron orbit to a drift orbit. As 
described above, the direction of rotation of the drift orbit is determined by a competition 
between the ExB drift and the gradient drift. If the ExB drift dominates, the drift orbit 
rotates in the diamagnetic direction. If the gradient drift dominates, the drift orbit rotates in 

20 the counterdiamagnetic direction. This is shown in Figs. 18A and 18B. Fig. 18A shows a 
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transition from a betatron orbit to a drift orbit due to a 180° collision, which occurs at the 
point 172. The drift orbit continues to rotate in the diamagnetic direction because the ExB 
drift dominates. Fig. 18B shows another 180° collision, but in this case the electrostatic field 
is weak and the gradient drift dominates. The drift orbit thus rotates in the 
counterdiamagnetic direction. 

The direction of rotation of the drift orbit determines whether it is confined or not. A 
particle moving in a drift orbit will also have a velocity parallel to the FRC axis. The time it 
takes the particle to go from one end of the FRC to the other, as a result of its parallel motion, 
is called transit time; thus, the drift orbits reach an end of the FRC in a time of the order of 
the transit time. As shown in connection with Fig. 15 A, the Lorentz force at the ends is 
confining only for drift orbits rotating in the diamagnetic direction. After a transit time, 
therefore, ions in drift orbits rotating in the counterdiamagnetic direction are lost. 

This phenomenon accounts for a loss mechanism for ions, which is expected to have 
existed in all FRC experiments. In fact, in these experiments, the ions carried half of the 
current and the electrons carried the other half. In these conditions the electric field inside 
the plasma was negligible, and the gradient drift always dominated the ExB drift. Hence, 
all the drift orbits produced by large-angle collisions were lost after a transit time. These 
experiments reported ion diffusion rates that were faster than those predicted by classical 
diffusion estimates. 
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If there is a strong electrostatic field, the ExB drift dominates the gradient drift, and 
the drift orbits rotate in the diamagnetic direction. This was shown above in connection with 
Fig. 18 A. When these orbits reach the ends of the FRC, they are reflected back into the 
region of closed field lines by the Lorentz force; thus, they remain confined in the system. 

The electrostatic fields in the colliding beam system may be strong enough, so that the 

ExB drift dominates the gradient drift. Thus, the electrostatic field of the system would 
avoid ion transport by eliminating this ion loss mechanism, which is similar to a loss cone in 
a mirror device. 

Another aspect of ion diffusion can be appreciated by considering the effect of small- 
angle, electron-ion collisions on betatron orbits. Fig. 19A shows a betatron orbit 112; Fig. 
19B shows the same orbit 112 when small-angle electron-ion collisions are considered 174; 
Fig. 19C shows the orbit of Fig. 19B followed for a time that is longer by a factor of ten 176; 
and Fig. 19D shows the orbit of Fig. 19B followed for a time longer by a factor of twenty 
178. It can be seen that the topology of betatron orbits does not change due to small-angle, 
electron-ion collisions; however, the amplitude of their radial oscillations grows with time. 
In fact, the orbits shown in Figs. 19A to 19D fatten out with time, which indicates classical 
diffusion. 
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Theory of the Invention 

For the purpose of modeling the invention, a one-dimensional equilibrium model for 
the colliding beam system is used, as shown in Fig. 10. The results described above were 
drawn from this model. This model shows how to derive equilibrium expressions for the 
particle densities, the magnetic field, the electric field, and the electric potential. The 
equilibrium model presented herein is valid for a plasma fuel with one type of ions (e.g., in a 
D-D reaction) or multiple types of ions (e.g., D-T, D-He 3 , and /?-B H ). 

Vlasov-Maxwell Equations 

Equilibrium solutions for the particle density and the electromagnetic fields in a FRC 
are obtained by solving self-consistently the Vlasov-Maxwell equations: 


9fj 


+ (v'V)fj+— E +—xB V v /. =0 


(5) 


dt 


m , c 



(6) 



(7) 



(8) 


j 


VB = 0, 


(9) 
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where j = e, i and i = 1, 2,. . . for electrons and each species of ions. In equilibrium, all 
physical quantities are independent of time (i.e., d/dt = 0). To solve the Vlasov-Maxwell 
equations, the following assumptions and approximations are made: 

(a) All the equilibrium properties are independent of axial position z (i.e., dldz = 0). 
This corresponds to considering a plasma with an infinite extension in the axial direction; 
thus, the model is valid only for the central part 88 of a FRC. 

(b) The system has cylindrical symmetry. Hence, all equilibrium properties do not 
depend on 0 (i.e., d/dQ = 0). 

(c) The Gauss law, Eq. 8, is replaced with the quasi-neutrality condition: Y, j n j e j = 0 . 

By assuming infinite axial extent of the FRC and cylindrical symmetry, all the equilibrium 
properties will depend only on the radial coordinate r. For this reason, the equilibrium model 
discussed herein is called one-dimensional. With these assumptions and approximations, the 
Vlasov-Maxwell equations reduce to: 


(v.V)/. + ^.V v / ; .+-^[vx4v v / ; =0 
m. mfi 


(10) 



j 


(ii) 


(12) 


a 
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Rigid Rotor Distributions 

To solve Eqs. 10 through 12, distribution functions must be chosen that adequately 
describe the rotating beams of electrons and ions in a FRC. A reasonable choice for this 
purpose are the so-called rigid rotor distributions, which are Maxwellian distributions in a 
uniformly rotating frame of reference. Rigid rotor distributions are functions of the constants 
of motion: 


fj(r,v) = 


k 2kT u 


«,(0)exp 


£ j-<»j p j 


(13) 


where ntj is particle mass, v is velocity, 7} is temperature, n/0) is density at r = 0, and oj is a 
constant. The constants of the motion are 

m i 2 

Sj = —*-v + e v O (for energy) and 


Pj = ntj (xv y - yv x )+ — *F (for canonical angular momentum), 
c 


where O is the electrostatic potential and is the flux function. The electromagnetic fields 
are 

ao 

E r = (electric field) and 

dr 

B z = (magnetic field). 
r dr 

Substituting the expressions for energy and canonical angular momentum into Eq. 13 yields 
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(14) 


where 


| v - 3 j x r | 2 = (v, + ya> j } + (v y - *<y y ) 2 + v) 


and 


»>('') = » y (0)exp 


7 . 


<D ^ 


J 


v 


^<y 2 r 2 


(15) 


That the mean velocity in Eq. 14 is a uniformly rotating vector gives rise to the name rigid 
rotor. One of skill in the art can appreciate that the choice of rigid rotor distributions for 
describing electrons and ions in a FRC is justified because the only solutions that satisfy 
Vlasov's equation (Eq. 10) are rigid rotor distributions (e.g., Eq. 14). A proof of this 
assertion follows: 
Proof 

We require that the solution of Vlasov's equation (Eq. 10) be in the form of a drifted 
Maxwellian: 


w,(r)exp 


27. (r) 


(16) 


i.e., a Maxwellian with particle density «,(r), temperature Tj(r) , and mean velocity Uj(r) 
that are arbitrary functions of position. Substituting Eq. 16 into the Vlasov's equation (Eq. 
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10) shows that (a) the temperatures Tj(r) must be constants; (b) the mean velocities Uj{r) 
must be uniformly rotating vectors; and (c) the particle densities iij(r) must be of the form of 
Eq. 15. Substituting Eq. 16 into Eq. 10 yields a third-order polynomial equation in v : 

/ \ m.[v -u .) . v m Av — u ;Y , x 
v-v(lnn,)+ A t W ; (v-V)u y+ A 2t2 <> (v -Vfe... 

... + ^£.(v--^^[?x5].(v-« y ) = 0. 
Grouping terms of like order in v yields 



For this polynomial equation to hold for all v , the coefficient of each power of v must 
vanish. 

The third-order equation yields Tj{r) = constant. 
The second-order equation gives 
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(du, du v du,i v A 


v-Vuj-v = (v x v y v z ) 


dx dx dx 

du x du y du 2 

dy dy dy 

du x du y du z 

< dz dz dz 


2 dU x l^V 2^2 

■ v] — - + v, — - + v] — - + v v 
* dx y dy 1 dz x y 


... + V V 


(du z 


) + v. 


du y } 

^ dx 


[dy 

+ — - 

dz j 


— - + — - 
{dx dy ) 


= 0. 


For this to hold for all v , we must satisfy 


^ = ^ = ^ = 0 and 

dx dy dz 

which is solved generally by 


(du v duA (du, duA (du. du 


dx dy 


K dx dz 


dy dz j 


= 0, 


(17) 


a j(f) = (ajxr)+u 0j 

In cylindrical coordinates, take u 0J = 0 and <y y = o)jZ, which corresponds to injection 
perpendicular to a magnetic field in the z direction. Then, u } (r) = q)j rO . 

The zero order equation indicates that the electric field must be in the radial direction, 
i.e., E = E r r. 

The first-order equation is now given by 

y.y(l nnj y^. V u j -u J )-^v-E + ^{vxB)u j =0. (18) 
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The second term in Eq. 18 can be rewritten with 


du r 

du 0 

du z > 

V 

dr 

dr 

dr ; 


1 du r 

1 du g 

1 du 2 

u e 

r dO 

r 80 

r 89 


du r 

du e 

8u 2 


dz 

dz 

dz j 

V J 


r 0 

0) 

j 

°1 

f 0 ^ 



0 

0 

0 

°>J r 


(19) 

,0 

0 


I o J 




The fourth term in Eq. 18 can be rewritten with 

(vxB}u j =v\Bxu J )=v\(yxA)xu J )=v^ (^(rA^x^cD.re) 

\\r dr ) 


= v-cOjj-(rA 0 )r 


Using Eqs. 19 and 20, the first-order Eq. 18 becomes 


— (inn.) Ltfr-^-E r +^^ — ( rA e ( r)) = 0 


8r 


cTj 8r 


The solution of this equation is 

rtj(r)= «,.(o)exp 

where E r = -dO/dr and «j(0) is given by 

",(°)="yo ex P 
Here, /ijo is the peak density at r 0 . 


mjQ)jr 2 e,<X>(r) ejO)jrA 0 (r) 


2Tj 


H h - 


27\ 


(20) 


(21) 


(22) 
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Solution of Vlasov-Maxwell Equations 

Now that it has been proved that it is appropriate to describe ions and electrons by 
rigid rotor distributions, the Vlasov's equation (Eq. 10) is replaced by its first-order moments, 


i.e., 


n m rco = n e ■ 

j j j j j 


rco, 
E + — J -B. 


dn, 
-T 1 1 


J dr ' 


which are conservation of momentum equations. The system of equations to obtain 
equilibrium solutions reduces to: 


rtjnijrcOj = n y e 7 - 


rco, 
E. + — J -B. 


dn., 


1 dr 


j = e,i = \,2,- 


d 1 5¥ dB 2 An . An: ^ 

= = J e = r Z_, n i e i CO i 

dr r dr dr c c , 


2>,e,=0. 


(23) 


(24) 


(25) 


(26) 


Solution for Plasma with One Type of Ion 

Consider first the case of one type of ion fully stripped. The electric charges are given 
by ej = -e.Ze. Solving Eq. 24 for E r with the electron equation yields 


m 2 rco. 


■ ~„ „ T. dn, 

E =—rco l -B 

e c en e dr 


and eliminating E r from the ion equation yields 

1 rflogK , _ Z,e {co i -G> e ) B _ ^Je_ ]_ d\ogn e + w,-6>, 2 + mZiCo] 


r dr 


Z r dr 


T, 


(27) 


(28) 
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Differentiating Eq. 28 with respect to r and substituting Eq. 25 for dB z /dr yields 

dB, An 


dr 


■n e er(a> t -G> e ) and Z,n, =n e , 


with T e = Tj = constant, and a*, co e , constants, obtaining 

1 d 1 </log«, 4;m e Z,e 2 {co l -Q) e ) 2 Z t T e 1 d 1 d\ogn e 


r dr r dr 


T, 


T r dr r dr 


(29) 


The new variable £ is introduced: 


\ d \ d Id'' 


2r 


r dr r dr r 0 d % 


(30) 


Eq. 29 can be expressed in terms of the new variable 


d 2 log«,. = 4^n e Z,e 2 r 0 4 (a?. -0) e f Z,T e </ 2 log» e 
rf 2 £ 7) c 2 r, rf 2 £ 


(31) 


Using the quasi-neutrality condition, 


10 


n e = Z,n t 


d 2 logw e _ d 2 log», 
</ 2 £ ~~d 2 £~ 


yields 


</ 2 ^ (?;+z,r,) c 2 


= 


l 


4aZ?e 2 (<y,-<y e ) 2 


*i0 


4^ e0 e 2 (<w,.-<a e ) 2 


-8 


>b 1 n 


-M —.(32) 


K&rJ n iQ 


Here is defined 
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r n Ar ■ 


T. 

T e+ — 

Z. 


4mt e0 e 


\0). 


(33) 


where the meaning of Ar will become apparent soon. If Ni = njnuk where w,o is the peak 
density at r = ro, Eq. 32 becomes 


(34) 


Using another new variable, 


J - , yields — + 
Ar d x 


= -2N n 


the solution to which is 


N; = 


1 


cosh 2 [x -ZoY 

where xo - X( r o) because of the physical requirement that Af,-(r 0 ) = 1. 
Finally, the ion density is given by 


n 


n. = 


iO 


•i0 


cosh 2 2 


I 2j 


cosh' 


(35) 


The significance of r 0 is that it is the location of peak density. Note that n i (0) = n i (V2r 0 ). 

With the ion density known, B z can be calculated using Eq. 11, and E r can be calculated using 
Eq. 27. 
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The electric and magnetic potentials are 

Q> = -£ =r o E r {r')dr' and 

A g = - [~ r r'B z (r')Jr* ¥ = rA g (flux function) (36) 

y «r'=0 

Taking r = V2r 0 to be the radius at the wall (a choice that will become evident when the 

expression for the electric potential O(r) is derived, showing that at r = 4lr Q the potential is 
zero, i.e., a conducting wall at ground), the line density is 


N < = Z > N ' = L — rs_„2\ = 2m ^ r tanh ± 

cosh 


. . . = 2nn e0 r 0 Ar (because r 0 » Ar ) (37) 
Thus, Ar represents an "effective thickness." In other words, for the purpose of line density, 
the plasma can be thought of as concentrated at the null circle in a ring of thickness Ar with 
constant density n e0 . 

The magnetic field is 

B z {r) = BM~~ C/r'nM^ ~<o e ). (38) 
The current due to the ion and electron beams is 

1 e ~ J, Je dr ~ — Je = n 0 er[p s -(O e ). 

Using Eq. 39, the magnetic field can be written as 


(39) 
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B 2 {r) = B,(0)-—I e -—I e tW-—^ = -B n ~—I a tanh ^ 


In w 2n 


r 0 Ar 


'0 *8 

c 


r 0 Ar 


(40) 


In Eq. 40, 


B z (0) = -B 0 +—I e and 
c 


B 2 (j2r o )=-B o -^I 0 


10 


If the plasma current Ig vanishes, the magnetic field is constant, as expected. 

These relations are illustrated in Figs. 20A through 20C. Fig. 20A shows the external 
magnetic field B Q 180. Fig. 20B shows the magnetic field due to the ring of current 182, the 

magnetic field having a magnitude of {2%lc)Ie. Fig. 20C shows field reversal 184 due to the 
overlapping of the two magnetic fields 180,182. 
The magnetic field is 


B 2 {r) = ~B 0 


. 2x1,, , r 2 -r 2 
1 + ^-tanh- 


cB n 


r 0 Ar 


= -B, 


l + V/?tanh 


f 2 2 


(41) 


using the following definition for /?: 


2n I 6 _ A^efc?, -co e ) _ 2k r^e{co i - co e ) 


c B n 


cB n 


= — n 


e0 


B n 


^24i 

C 


4mt e0 e 2 


2 c» e0 e(a>,-a> e ) 


CO; - CO. 


B n 
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10 


Bl 


(42) 


With an expression for the magnetic field, the electric potential and the magnetic flux 
can be calculated. From Eq. 27, 


rco T e d\nn m 2 d<S> 
E. = -B, — — - + — rco; = 


(43) 


c 2 e dr e dr 
Integrating both sides of Eq. 28 with respect to r and using the definitions of electric 
potential and flux function, 

O^-PX^' and (~ r BAr % )r % dr\ (44) 

which yields 


e e « e (0) e 2 


(45) 


Now, the magnetic flux can be calculated directly from the expression of the magnetic 
field (Eq.41): 


r- r 

J-'=0 


2 2 

1 + V/Ttanh °- r'dr' 


r 0 Ar 


B 0 r 2 B o j0 


r 0 Ar 



( 2 1 ~\ 


( r > 


log 

cosh r r ° 

-log 

cosh -2- 


I ro& ) 

l ArJ 



2 0 4 % e fo; 


(46) 


15 Substituting Eq. 46 into Eq. 45 yields 
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(47) 


Using the definition of /?, 


C "* 


CO 

c 


... = 4- 


®e ( n eO T e+ n iO T i) 


0) —CO 
i e 


(48) 


Finally, using Eq. 48, the expressions for the electric potential and the flux function become 


*F(r) = - 


By 


CO; - a. 


r n e0 T e +n i0 T^ 


ln-^ and 
».(0) 


<D(r) = 


e 


CO - - co 

i e 


In 

fl c (0) c 2 e c 


(49) 


(50) 


Relationship Between a* and <u 


An expression for the electron angular velocity co e can also be derived from Eqs. 24 
through 26. It is assumed that ions have an average energy J4/w,{rfi*) 2 , which is determined by 
the method of formation of the FRC. Therefore, a* is determined by the FRC formation 
method, and co e can be determined by Eq. 24 by combining the equations for electrons and 
ions to eliminate the electric field: 


T 2 2 1 n ? er ( \ n ^ an, m an. 

c dr dr 


(51) 


15 Eq. 25 can then be used to eliminate (^ - co e ) to obtain 
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[n e mrco] + n^rcof ] = + £ n j T j 


dr 

Eq. 52 can be integrated from r = 0 to r B = V2r 0 . Assuming ro/Ar » 1, the density is very 
small at both boundaries arid B 2 = -B 0 (l ± Carrying out the integration shows 


(52) 


[n e0 mo>l + « 1 . 0 i Wf o I 2 ]r 0 Ar = -±[%n{n e0 T e + n^JJi . (53) 


Using Eq. 33 for Ar yields an equation for co e \ 


^ 2 =O 0 fo-<O> (54) 


where Q 0 = - 


Some limiting cases derived from Eq. 54 are: 


1. a* = 0 and co e = — 


'o . 


mc 


2. a> e = 0 and co i = Q 0 ; and 


Zm 2 2 j 

3. co e « co i and co e ~ co i 

m, 


In the first case, the current is carried entirely by electrons moving in their 
diamagnetic direction (co e < 0). The electrons are confined magnetically, and the ions are 
confined electrostatically by 


DOCSOCl: 147737.1 
703538^039 K2R 


44 


PATENT 
703538.4039 

T. dn. < 0 for r>r 0 

£=— 1 •- °. (55) 

Zen, dr >0forr<r 0 

In the second case, the current is carried entirely by ions moving in their diamagnetic 
direction (an > 0). If <uj is specified from the ion energy Vmfrcoi) 2 , determined in the 
formation process, then co e - 0 and Q 0 = co i identifies the value of B 0 , the externally applied 

magnetic field. The ions are magnetically confined, and electrons are electrostatically 
confined by 

T e dn e >0forr>r o 

L r J" . A ~ ^ • (56) 

en e dr < 0 for r < r 0 

In the third case, co e > 0 and Q 0 > co x . Electrons move in their counter diamagnetic 

direction and reduce the current density. From Eq. 33, the width of the distribution n t {r) is 
increased; however, the total current/unit length is 

*' = llo J ° dr = ^t 6 ^ ~ WherC (5?) 

N e = [[2nrdrn e = 2xr 0 Arn e0 . (58) 

Here, r B = 4lr Q and r 0 Ar oc {co i -o) e )" 1 according to Eq. 33. The electron angular 
velocity co e can be increased by tuning the applied magnetic field 2? 0 . This does not change 
either 1$ or the maximum magnetic field produced by the plasma current, which is 
B o VP = ( 2n /c Ve ■ However, it does change Ar and, significantly, the potential O. The 
maximum value of O is increased, as is the electric field that confines the electrons. 
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Tuning the Magnetic Field 

In Figs. 21 A-D, the quantities njn& 186, B z /(b 0 Vp) 188, a>AD 0 190, and 192 
are plotted against r/r 0 194 for various values of j5 0 . The values of potential and flux are 
normalized to O 0 = 20(7; + Ttye and To = (c/co,)0 0 . A deuterium plasma is assumed with the 
following data: n# = n i0 = 10 15 cm" 3 ; r 0 = 40 cm; ^{rofi*) 2 = 300 keV; and T e = Ti= 100 
keV. For each of the cases illustrated in Fig. 21, co i =1.35 x 10 7 s" 1 , and co e is determined 
from Eq. 54 for various values of Bq: 


Plot 

applied magnetic field (Bo) 

electron angular velocity (a> e ) 

154 

B 0 = 2.77 kG 

o) e = 0 

156 

£ 0 = 5.15kG 

co e = 0.625 x 10V 

158 

B 0 = 15.5 kG 

co e =\.U x 10 7 s" 1 


The case of co e = -0% and Bo = 1.385 kG involves magnetic confinement of both 
electrons and ions. The potential reduces to O/O 0 = m i {ro^) 2 l[? i 0{T e + 7})], which is 
negligible compared to the case co e = 0. The width of the density distribution Ar is reduced 
by a factor of 2, and the maximum magnetic field B 0 ^[fi is the same as for co e = 0. 


Solution for Plasmas of Multiple Types of Ions 

This analysis can be carried out to include plasmas comprising multiple types of ions. 
Fusion fuels of interest involve two different kinds of ions, e.g., D-T, D-He 3 , and H-B 11 . 
The equilibrium equations (Eqs. 24 through 26) apply, except that j = e, 1,2 denotes 
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electrons and two types of ions where Z\ = 1 in each case and Zi = Z = 1, 2, 5 for the above 
fuels. The equations for electrons and two types of ions cannot be solved exactly in terms of 
elementary functions. Accordingly, an iterative method has been developed that begins with 
an approximate solution. 

The ions are assumed to have the same values of temperature and mean velocity V[ = 
re®. Ion-ion collisions drive the distributions toward this state, and the momentum transfer 
time for the ion-ion collisions is shorter than for ion-electron collisions by a factor of an order 
of 1000. By using an approximation, the problem with two types of ions can be reduced to a 
single ion problem. The momentum conservation equations for ions are 


-n x m x rco x - n x e 


c 


-T x *± and 
1 dr 


n 2 m 2 rcol = n 2 Ze 


E r + — 2 -B, 


rr, dn 7 


dr 


(59) 


(60) 


In the present case, T\ = T 2 and co\ = o^. Adding these two equations results in 


W =n,(Z)e 


-T, 


dm 


1 dr' 


(61) 


where n\ = n\ + n 2 ; gh = 0)i = coi, T t = T\- Ti\ "/(/«,.) =n\ni\ + nitni, and H;(Z) ="i + n 2 Z. 
The approximation is to assume that (/«, } and (z) are constants obtained by 

replacing «i(r) and ni(r) by «io and «2o, the maximum values of the respective functions. The 
solution of this problem is now the same as the previous solution for the single ion type, 
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except that (Z) replaces Z and ) replaces m,. The values of n\ and «2 can be obtained 

from n\ + Ai2 = 1/ and /ij + Zn 2 = " e = • I* can ^ e appreciated that «i and ni have the 
same functional form. 

Now the correct solution can be obtained by iterating the equations: 

i^^ Ckx ^^m.^i^^W^ md (62) 

d$ 1 0 1 7; 5 0 t; </£ 7;. v ' 

m ^ feziO Ml _ S + afeaZ , (63) 

where 


j^ilW, Af 2 =^, £=^, Q,=^,. and Q 2 =^. 

«10 W 20 2r 0 W 1 C ™2 C 

The first iteration can be obtained by substituting the approximate values of B z (^) and 
N e (Q in the right hand sides of Eqs. 62 and 63 and integrating to obtain the corrected values 
ofwi(r), n 2 (r), and£ z (r). 

Calculations have been carried out for the data shown in Table 1, below. Numerical 
results for fusion fuels are shown in Figs. 22 A-D through 24 A-D wherein the quantities 
ni/nio 206, OAD 0 208, and 47% 210 are plotted against r/r 0 204. Figs. 22 A-D shows the 
first approximation (solid lines) and the final results (dotted lines) of the iterations for D-T 
for the normalized density of D 196, the normalized density of T 198, the normalized electric 
potential 200, and the normalized flux 202. Figs. 23 A-D show the same iterations for D-He 3 
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for the normalized density of D 212, the normalized density of He 3 214, the normalized 
electric potential 216, and the normalized flux 218. Figs. 24 A-D show the same iterations 
for p-B n for the normalized density ofp 220, the normalized density of B 11 222, the 
normalized electric potential 224, and the normalized flux 226. Convergence of the iteration 
5 is most rapid for D-T. In all cases the first approximation is close to the final result. 


Table 1: Numerical data for equilibrium calculations for different fusion fuels 


Quantity 

Units 

D-T 

D-He 3 

p-B n 


cm" 3 

10 15 

10 15 

10 15 

"10 

cm" 3 

0.5 x 10 15 

1/3 x 10 15 

0.5 x 10 15 

"20 

cm' 3 

0.5 x 10 15 

1/3 x 10 15 

10 14 

Vi = V 2 

cm 
sec 

0.54 x 10 9 

0.661 x 10 9 

0.764 x 10 9 

1 2 

—m,v, 

2 1 1 

keV 

300 

450 

300 

1 2 

keV 

450 

675 

3300 

0\= CO\ = (Ol 

rad/s 

1.35 x 10 7 

1.65 x 10 7 

1.91 x 10 7 

r 0 

cm 

40 

40 

40 

Bo 

kG 

5.88 

8.25 

15.3 

(Zi) 

None 

1 

3/2 

1.67 


m p 

5/2 

5/2 

2.67 
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V^uantity 

units 

LJ-1 

D-rle 

- nil 

n 0 = {Zi)eB ° 
\ m i) c 

rad/s 

o o c . . 1 /\7 

2.35 x 10 

4.95 x 10 

9.55 x 10 7 

CO- 
CO. = CO; 1 

rad/s 

0 x 10 7 

v. J / J A 1 \J 

1 1 x 10 7 

1 s? * 1 n 7 


keV 

96 

170 

82 

Ti 

keV 

100 

217 

235 

roAr 

cm 2 

114 

203 

313 

P 

None 

228 

187 

38.3 


Structure of the Containment System 


Fig. 25 illustrates a preferred embodiment of a containment system 300 according to 
the present invention. The containment system 300 comprises a chamber wall 305 that 
defines therein a confining chamber 310. Preferably, the chamber 310 is cylindrical in shape, 
with principle axis 315 along the center of the chamber 310. For application of this 
containment system 300 to a fusion reactor, it is necessary to create a vacuum or near vacuum 
inside the chamber 310. Concentric with the principle axis 315 is a betatron flux coil 320, 
located within the chamber 310. The betatron flux coil 320 comprises an electrical current 
carrying medium adapted to direct current around a long coil, as shown, which preferably 
comprises parallel winding multiple separate coils, and most preferably parallel windings of 
about four separate coils, to form a long coil. Persons skilled in the art will appreciate that 
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current through the betatron coil 320 will result in a magnetic field inside the betatron coil 
320, substantially in the direction of the principle axis 315. 

Around the outside of the chamber wall 305 is an outer coil 325. The outer coil 325 
produce a relatively constant magnetic field having flux substantially parallel with principle 
5 axis 315. This magnetic field is azimuthally symmetrical. The approximation that the 

magnetic field due to the outer coil 325 is constant and parallel to axis 315 is most valid away 
from the ends of the chamber 310. At each end of the chamber 3 10 is a mirror coil 330. The 
mirror coils 330 are adapted to produce an increased magnetic field inside the chamber 310 at 
each end, thus bending the magnetic field lines inward at each end. (See Figs. 8 and 10.) As 

10 explained, this bending inward of the field lines helps to contain the plasma 335 in a 

containment region within the chamber 310 generally between the mirror coils 330 by 
pushing it away from the ends where it can escape the containment system 300. The mirror 
coils 330 can be adapted to produce an increased magnetic field at the ends by a variety of 
methods known in the art, including increasing the number of windings in the mirror coils 

15 330, increasing the current through the mirror coils 330, or overlapping the mirror coils 330 

with the outer coil 325. 

The outer coil 325 and mirror coils 330 are shown in Fig. 25 implemented outside the 
chamber wall 305; however, they may be inside the chamber 310. In cases where the 
chamber wall 305 is constructed of a conductive material such as metal, it may be 

20 advantageous to place the coils 325, 330 inside the chamber wall 305 because the time that it 
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takes for the magnetic field to diffuse through the wall 305 may be relatively large and thus 
cause the system 300 to react sluggishly. Similarly, the chamber 310 may be of the shape of 
a hollow cylinder, the chamber wall 305 forming a long, annular ring. In such a case, the 
betatron flux coil 320 could be implemented outside of the chamber wall 305 in the center of 
that annular ring. Preferably, the inner wall forming the center of the annular ring may 
comprise a non-conducting material such as glass. As will become apparent, the chamber 
310 must be of sufficient size and shape to allow the circulating plasma beam or layer 335 to 
rotate around the principle axis 315 at a given radius. 

The chamber wall 305 may be formed of a material having a high magnetic 
permeability, such as steel. In such a case, the chamber wall 305, due to induced 
countercurrents in the material, helps to keep the magnetic flux from escaping the chamber 
310, "compressing" it. If the chamber wall were to be made of a material having low 
magnetic permeability, such as plexiglass, another device for containing the magnetic flux 
would be necessary. In such a case, a series of closed-loop, flat metal rings could be 
provided. These rings, known in the art as flux delimiters, would be provided within the 
outer coils 325 but outside the circulating plasma beam 335. Further, these flux delimiters 
could be passive or active, wherein the active flux delimiters would be driven with a 
predetermined current to greater facilitate the containment of magnetic flux within the 
chamber 310. Alternatively, the outer coils 325 themselves could serve as flux delimiters. 
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As explained above, a circulating plasma beam 335, comprising charged particles, 
may be contained within the chamber 310 by the Lorentz force caused by the magnetic field 
due to the outer coil 325. As such, the ions in the plasma beam 335 are magnetically 
contained in large betatron orbits about the flux lines from the outer coil 325, which are 
parallel to the principle axis 315. One or more beam injection ports 340 are also provided for 
adding plasma ions to the circulating plasma beam 335 in the chamber 310. In a preferred 
embodiment, the injector ports 340 are adapted to inject an ion beam at about the same radial 
position from the principle axis 315 where the circulating plasma beam 335 is contained (i.e., 
around the null surface). Further, the injector ports 340 are adapted to inject ion beams 350 
(See Fig. 28) tangent to and in the direction of the betatron orbit of the contained plasma 
beam 335. 

Also provided are one or more background plasma sources 345 for injecting a cloud 
of non-energetic plasma into the chamber 310. In a preferred embodiment, the background 
plasma sources 345 are adapted to direct plasma 335 toward the axial center of the chamber 
310. It has been found that directing the plasma this way helps to better contain the plasma 
335 and leads to a higher density of plasma 335 in the containment region within the chamber 
310. 
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Formation of the FRC 

Conventional procedures used to form a FRC primarily employ the theta pinch- field 
reversal procedure. In this conventional method, a bias magnetic field is applied by external 
coils surrounding a neutral gas back-filled chamber. Once this has occurred, the gas is 
ionized and the bias magnetic field is frozen in the plasma. Next, the current in the external 
coils is rapidly reversed and the oppositely oriented magnetic field lines connect with the 
previously frozen lines to form the closed topology of the FRC (see Fig. 8). This formation 
process is largely empirical and there exists almost no means of controlling the formation of 
the FRC. The method has poor reproducibility and no tuning capability as a result. 

In contrast, the FRC formation methods of the present invention allow for ample 
control and provide a much more transparent and reproducible process. In fact, the FRC 
formed by the methods of the present invention can be tuned and its shape as well as other 
properties can be directly influenced by manipulation of the magnetic field applied by the 
outer field coils 325. Formation of the FRC by methods of the present inventions also results 
in the formation of the electric field and potential well in the manner described in detail 
above. Moreover, the present methods can be easily extended to accelerate the FRC to reactor 
level parameters and high-energy fuel currents, and advantageously enables the classical 
confinement of the ions. Furthermore, the technique can be employed in a compact device 
and is very robust as well as easy to implement - all highly desirable characteristics for 
reactor systems. 
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In the present methods, FRC formation relates to the circulating plasma beam 335. It 
can be appreciated that the circulating plasma beam 335, because it is a current, creates a 
poloidal magnetic field, as would an electrical current in a circular wire. Inside the 
circulating plasma beam 335, the magnetic self-field that it induces opposes the externally 
applied magnetic field due to the outer coil 325. Outside the plasma beam 335, the magnetic 
self-field is in. the same direction as the applied magnetic field. When the plasma ion current 
is sufficiently large, the self-field overcomes the applied field, and the magnetic field reverses 
inside the circulating plasma beam 335, thereby forming the FRC topology as shown in Figs. 
8 and 10. 

The requirements for field reversal can be estimated with a simple model. Consider 
an electric current I P carried by a ring of major radius r 0 and minor radius a « r 0 . The 
magnetic field at the center of the ring normal to the ring is B p = 2nI P /(cr 0 ). Assume that the 
ring current I P = N p e(Q^2n) is carried by N p ions that have an angular velocity fi 0 - For a 
single ion circulating at radius r 0 = Fb/Qo, &o = eitym/C is the cyclotron frequency for an 
external magnetic field B 0 . Assume V 0 is the average velocity of the beam ions. Field 
reversal is defined as 

N D eCl Q 

B p =^^>2B 09 (64) 
r 0 c 


which implies that N p >2 ro/ar,, and 


p 7ta : 
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where a t = e 2 l mp 2 = 1.57 x 10" 16 cm and the ion beam energy is \myl . In the one- 
dimensional model, the magnetic field from the plasma current is B p = (2n/c)i py where i p is 
current per unit of length. The field reversal requirement is i p > eVo/nr 0 ai = 0.225 kA/cm, 
where Bo = 69.3 G and \m^l =100 eV. For a model with periodic rings and B z is averaged 
over the axial coordinate (B 2 ) = {l7tlc){l p Is) (s is the ring spacing), if s = r 0 , this model 

would have the same average magnetic field as the one dimensional model with i p = Ipls. 
Combined Beam/Betatron Formation Technique 

A preferred method of forming a FRC within the confinement system 300 described 
above is herein termed the combined beam/betatron technique. This approach combines low 
energy beams of plasma ions with betatron acceleration using the betatron flux coil 320. 

The first step in this method is to inject a substantially annular cloud layer of 
background plasma in the chamber 310 using the background plasma sources 345. Outer coil 
325 produces a magnetic field inside the chamber 310, which magnetizes the background 
plasma. At short intervals, low energy ion beams are injected into the chamber 310 through 
the injector ports 340 substantially transverse to the externally applied magnetic field within 
the chamber 310. As explained above, the ion beams are trapped within the chamber 310 in 
large betatron orbits by this magnetic field. The ion beams may be generated by an ion 
accelerator, such as an accelerator comprising an ion diode and a Marx generator, (see R.B. 
Miller, An Introduction to the Physics of Intense Charged Particle Beams, (1982)). As one of 
skill in the art can appreciate, the externally applied magnetic field will exert a Lorentz force 
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on the injected ion beam as soon as it enters the chamber 310; however, it is desired that the 
beam not deflect, and thus not enter a betatron orbit, until the ion beam reaches the 
circulating plasma beam 335. To solve this problem, the ion beams are neutralized with 
electrons and directed through a substantially constant unidirectional magnetic field before 
5 entering the chamber 310. As illustrated in Fig. 26, when the ion beam 350 is directed 

through an appropriate magnetic field, the positively charged ions and negatively charged 
electrons separate. The ion beam 350 thus acquires an electric self-polarization due to the 
magnetic field. This magnetic field may be produced by, e.g., a permanent magnet or by an 
electromagnet along the path of the ion beam. When subsequently introduced into the 

10 confinement chamber 310, the resultant electric field balances the magnetic force on the 

beam particles, allowing the ion beam to drift undeflected. Fig. 27 shows a head-on view of 
the ion beam 350 as it contacts the plasma 335. As depicted, electrons from the plasma 335 
travel along magnetic field lines into or out of the beam 350, which thereby drains the beam's 
electric polarization. When the beam is no longer electrically polarized, the beam joins the 

1 5 circulating plasma beam 335 in a betatron orbit around the principle axis 3 1 5, as shown in 

Fig. 25. 

When the plasma beam 335 travels in its betatron orbit, the moving ions comprise a 
current, which in turn gives rise to a poloidal magnetic self-field. To produce the FRC 
topology within the chamber 310, it is necessary to increase the velocity of the plasma beam 
20 335, thus increasing the magnitude of the magnetic self-field that the plasma beam 335 
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causes. When the magnetic self- field is large enough, the direction of the magnetic field at 
radial distances from the axis 315 within the plasma beam 335 reverses, giving rise to a FRC. 
(See Figs. 8 and 10). It can be appreciated that, to maintain the radial distance of the 
circulating plasma beam 335 in the betatron orbit, it is necessary to increase the applied 
magnetic field from the outer coil 325 as the plasma beam 335 increases in velocity. A 
control system is thus provided for maintaining an appropriate applied magnetic field, 
dictated by the current through the outer coil 325. Alternatively, a second outer coil may be 
used to provide the additional applied magnetic field that is required to maintain the radius of 
the plasma beam's orbit as it is accelerated. 

To increase the velocity of the circulating plasma beam 335 in its orbit, the betatron 
flux coil 320 is provided. Referring to Fig. 28, it can be appreciated that increasing a current 
through the betatron flux coil 320, by Ampere's Law, induces an azimuthal electric field, E, 
inside the chamber 310. The positively charged ions in the plasma beam 335 are accelerated 
by this induced electric field, leading to field reversal as described above. When ion beams 
are added to the circulating plasma beam 335, as described above, the plasma beam 335 
depolarizes the ion beams. 

For field reversal, the circulating plasma beam 335 is preferably accelerated to a 
rotational energy of about 100 eV, and preferably in a range of about 75 eV to 125 eV. To 
reach fusion relevant conditions, the circulating plasma beam 335 is preferably accelerated to 
about 200 keV and preferably to a range of about 100 keV to 3.3 MeV. In developing the 
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necessary expressions for the betatron acceleration, the acceleration of single particles is first 
considered. The gyroradius of ions r = V/Qj will change because V increases and the applied 
magnetic field must change to maintain the radius of the plasma beam's orbit, ro = V/Q. c 

dV V dQ, 


9r = J_ 
dt~ Q 


dt Q ; dt 


= 0, 


(66) 


where 


dV_ = r^e_dB^ = eEg_^ e 1 

dt ntjC dt m ( m,c2m- 0 dt '. 


(67) 


and is the magnetic flux: 


V=[°B z 2m-dr = nr 0 2 {B 2 ), 


(68) 


where 



2 


( \ 



1- 


\ r o) 




(69) 


From Eq. 67, it follows that 


dt dt ' 

and (B z ) = -2B C + B 0 , assuming that the initial values of Br and B c are both Bq. Eq. 67 can 
be expressed as 


(70) 


dV_ 
dt 


2m : c u dt 


(71) 
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After integration from the initial to final states where \mVl = W Q and \mV 2 = W , the final 
values of the magnetic fields are: 


W 


B c =BJ— = 2.19kG 


(72) 


10 


15 


and 



+ 



-1 


= 10.7kG, 


(73) 


assuming B 0 = 69.3 G, 07 Wb = 1000, and r<Jr a = 2. This calculation applies to a collection of 
ions, provided that they are all located at nearly the same radius r 0 and the number of ions is 
insufficient to alter the magnetic fields. 

The modifications of the basic betatron equations to accommodate the present 
problem will be based on a one-dimensional equilibrium to describe the multi-ring plasma 
beam, assuming the rings have spread out along the field lines and the z-dependence can be 
neglected. The equilibrium is a self-consistent solution of the Vlasov-Maxwell equations that 
can be summarized as follows: 


(a) The density distribution is 


n = 


cosh' 


(74) 


which applies to the electrons and protons (assuming quasi neutrality); ro is the position of 
the density maximum; and Ar is the width of the distribution; and 
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(b) The magnetic field is 

27d D |V-r,M 

B 2 =-B c Manh — — £- , (75) 

c { r 0 Ar ) 

where B c is the external field produced by the outer coil 325. Initially, B c = B 0 : This solution 
satisfies the boundary conditions that r = r a and r = r b are conductors (5 n0 rmai = 0) and 
equipotentials with potential 0 = 0. The boundary conditions are satisfied if 
r 0 2 = (r fl 2 + rl )l 2 . r a = 10 cm and r 0 = 20 cm, so it follows that r b = 26.5 cm. I p is the 
plasma current per unit length. 

The average velocities of the beam particles are Vt = r 0 ^ and V e = r§co e , which are 
related by the equilibrium condition: 


where Q, = eBJ(mfi). Initially, it is assumed B c = B 0 ,(Oi = Q„ and co e = 0. (In the initial 
equilibrium there is an electric field such that the ExB and the VBxB drifts cancel. Other 
equilibria are possible according to the choice of B c .) The equilibrium equations are assumed 
to be valid if cch and B c are slowly varying functions of time, but r 0 = F}/Q,- remains constant. 
The condition for this is the same as Eq. 66. Eq. 67 is also similar, but the flux function *F 
has an additional term, i.e., ¥ = 7rr*(B z ) where 



(76) 


IbUb +— i p- 



a J 


(77) 


DOCSOC1: 147737.1 
703538^039 K2R 


61 


f 


10 


V 

2 


V 



1- 


\ r oJ 




and 

B z = -B F 

The magnetic energy per unit length due to the beam current is 


b.-bA 2 


&7I 


2 " " 


from which 


£ = r b 2 -r a 2 2n 2 r 2 


P r 2 + r 2 c 2 


(b 2 ) = b, + ^l p i p . 

The betatron condition of Eq. 70 is thus modified so that 


2- and 


^ = _ 2 ^- L f cdI " 


dt 


and Eq. 67 becomes: 


dV L _er k dB JL 


dt ni; c dt 


dt nr 2 dt 


e dB z e L p dl p 


Imp dt m i 2jir Q dt 


After integrating, 


Afl = -2B f 


r 2 -r 2 ' 
1+ " 2 ° 

1-51-1 

r o 
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(78) 


(79) 


(80) 


(81) 


(82) 


(83) 
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For W 0 = 100 eV and W= 100 keV, A5 Z = -7.49 kG . Integration of Eqs. 81 and 82 
determines the value of the magnetic field produced by the field coil: 

B e =B 0 i — =2.19kG 


(84) 


and 


(r ^ 

2 



AS,- 


So) 

{ >■: ) 


AB c = 25 kG . 


(85) 


If the final energy is 200 keV, B c = 3. 1 3 kG and B F = 34.5 kG. The magnetic energy 
B 2 

in the flux coil would be —Ttrjl = 172 kJ . The plasma current is initially 0.225 kA/cm 

corresponding to a magnetic field of 140 G, which increases to 10 kA/cm and a magnetic 
field of 6.26 kG. In the above calculations, the drag due to Coulomb collisions has been 
neglected. In the injection/trapping phase, it was equivalent to 0.38 volts/cm. It decreases as 
the electron temperature increases during acceleration. The inductive drag, which is 
included, is 4.7 volts/cm, assuming acceleration to 200 keV in 100 jas. 

The betatron flux coil 320 also balances the drag from collisions and inductance. The 
frictional and inductive drag can be described by the equation: 


dt 


i_ j_ 


_e L dl b 

m b 2nr 0 dt 


(86) 


where (7)//w,) < V b < (TJm). Here, V b is the beam velocity, T e and T t are electron and ion 
temperatures, h is the beam ion current, and 
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L = 0.01257r n 


In 


= 0.71/*H 


is the ring inductance. Also, r<> = 20 cm and a = 4 cm. 
The Coulomb drag is determined by 


3 [YfrnA 


ne 4 In A 


= 195/^sec 


4;zwe In A 


= 54.8// sec 


(87) 


To compensate the drag, the betatron flux coil 320 must provide an electric field of 1 .9 
volts/cm (0.38 volts/cm for the Coulomb drag and 1.56 volts/cm for the inductive drag). The 
magnetic field in the betatron flux coil 320 must increase by 78 Gausses to accomplish this, 
in which case Vt will be constant. The rise time of the current to 4.5 kA is 18 us, so that the 
magnetic field B F will increase by 1 .4 kG. The magnetic field energy required in the betatron 
flux coil 320 is 


B 2 

— x nr 2 F l = 394 Joules (/ = 1 1 5cm) . 

&7T 


(88) 


Betatron Formation Technique 

Another preferred method of forming a FRC within the confinement system 300 is 
herein termed the betatron formation technique. This technique is based on driving the 
betatron induced current directly to accelerate a circulating plasma beam 335 using the 
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betatron flux coil 320. A preferred embodiment of this technique uses the confinement 
system 300 depicted in Fig. 25, except that the injection of low energy ion beams is not 
necessary. 

As indicated, the main component in the betatron formation technique is the betatron 
flux coil 320 mounted in the center and along the axis of the chamber 310. Due to its 
separate parallel windings construction, the coil 320 exhibits very low inductance and, when 
coupled to an adequate power source, has a low LC time constant, which enables rapid ramp 
up of the current in the flux coil 320. 

Preferably, formation of the FRC commences by energizing the external field coils 
325, 330. This provides an axial guide field as well as radial magnetic field components near 
the ends to axially confine the plasma injected into the chamber 310. Once sufficient 
magnetic field is established, the background plasma sources 345 are energized from their 
own power supplies. Plasma emanating from the guns streams along the axial guide field and 
spreads slightly due to its temperature. As the plasma reaches the mid-plane of the chamber 
310, a continuous, axially extending, annular layer of cold, slowly moving plasma is 
established. 

At this point the betatron flux coil 320 is energized. The rapidly rising current in the 
coil 320 causes a fast changing axial flux in the coil's interior. By virtue of inductive effects 
this rapid increase in axial flux causes the generation of an azimuthal electric field E (see Fig. 
29), which permeates the space around the flux coil. By Maxwell's equations, this electric 
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field is directly proportional to the change in strength of the magnetic flux inside the coil, i.e.: 
a faster betatron coil current ramp-up will lead to a stronger electric field. 

The inductively created electric field couples to the charged particles in the plasma 
and causes a ponderomotive force, which accelerates the particles in the annular plasma layer. 
5 Electrons, by virtue of their smaller mass, are the first species to experience acceleration. 

The initial current formed by this process is, thus, primarily due to electrons. However, 
sufficient acceleration time (around hundreds of micro-seconds) will eventually also lead to 
ion current. Referring to Fig. 29, this electric field accelerates the electrons and ions in 
opposite directions. Once both species reach their terminal velocities, current is carried about 

1 0 equally by ions and electrons. 

As noted above, the current carried by the rotating plasma gives rise to a self magnetic 
field. The creation of the actual FRC topology sets in when the self magnetic field created 
by the current in the plasma layer becomes comparable to the applied magnetic field from the 
external field coils 325, 330. At this point magnetic reconnection occurs and the open field 

1 5 lines of the initial externally produced magnetic field begin to close and form the FRC flux 

surfaces (see Figs. 8 and 10). 

The base FRC established by this method exhibits modest magnetic field and particle 
energies that are typically not at reactor relevant operating parameters. However, the 
inductive electric acceleration field will persist, as long as the current in the betatron flux coil 

20 320 continues to increase at a rapid rate. The effect of this process is that the energy and total 
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magnetic field strength of the FRC continues to grow. The extent of this process is, thus, 
primarily limited by the flux coil power supply, as continued delivery of current requires a 
massive energy storage bank. However, it is, in principal, straightforward to accelerate the 
system to reactor relevant conditions. 
5 For field reversal, the circulating plasma beam 335 is preferably accelerated to a 

rotational energy of about 100 eV, and preferably in a range of about 75 eV to 125 eV. To 
reach fusion relevant conditions, the circulating plasma beam 335 is preferably accelerated to 
about 200 keV and preferably to a range of about 100 keV to 3.3 MeV. When ion beams are 
added to the circulating plasma beam 335, as described above, the plasma beam 335 

1 0 depolarizes the ion beams. 

Experiments - Beam Trapping and FRC Formation 

Experiment 1 : Propagating and trapping of a neutralized beam in a magnetic 
containment vessel to create an FRC. 

Beam propagation and trapping were successfully demonstrated at the following 

1 5 parameter levels : 

• Vacuum chamber dimensions: about 1 m diameter, 1.5 m length. 

• Betatron coil radius of 1 0 cm. 

• Plasma beam orbit radius of 20 cm. 
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• Mean kinetic energy of streaming beam plasma was measured to be about 100 eV, 
with a density of about 10 13 cm" 3 , kinetic temperature on the order of 10 eV and a 
pulse-length of about 20 jas. 

• Mean magnetic field produced in the trapping volume was around 100 Gauss, with a 
ramp-up period of 150 jis. Source: Outer coils and betatron coils. 

• Neutralizing background plasma (substantially Hydrogen gas) was characterized by a 
mean density of about 10 13 cm" 3 , kinetic temperature of less than 10 eV. 

The beam was generated in a deflagration type plasma gun. The plasma beam source 
was neutral Hydrogen gas, which was injected through the back of the gun through a special 
puff valve. Different geometrical designs of the electrode assembly were utilized in an 
overall cylindrical arrangement. The charging voltage was typically adjusted between 5 and 
7.5 kV. Peak breakdown currents in the guns exceeded 250,000 A. During part of the 
experimental runs, additional pre-ionized plasma was provided by means of an array of small 
peripheral cable guns feeding into the central gun electrode assembly before, during or after 
neutral gas injection. This provided for extended pulse lengths of above 25 jis. 

The emerging low energy neutralized beam was cooled by means of streaming 
through a drift tube of non-conducting material before entering the main vacuum chamber. 
The beam plasma was also pre-magnetized while streaming through this tube by means of 
permanent magnets. 
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The beam self-polarized while traveling through the drift tube and entering the 
chamber, causing the generation of a beam-internal electric field that offset the magnetic field 
forces on the beam. By virtue of this mechanism it was possible to propagate beams as 
characterized above through a region of magnetic field without deflection. 

Upon further penetration into the chamber, the beam reached the desired orbit 
location and encountered a layer of background plasma provided by an array of cable guns 
and other surface flashover sources. The proximity of sufficient electron density caused the 
beam to loose its self-polarization field and follow single particle like orbits, essentially 
trapping the beam. Faraday cup and B-dot probe measurements confirmed the trapping of the 
beam and its orbit. The beam was observed to have performed the desired circular orbit upon 
trapping. The beam plasma was followed along its orbit for close to 3 A of a turn. The 
measurements indicated that continued factional and inductive losses caused the beam 
particles to loose sufficient energy for them to curl inward from the desired orbit and hit the 
betatron coil surface at around the 3 A turn mark. To prevent this, the losses could be 
compensated by supplying additional energy to the orbiting beam by inductively driving the 
particles by means of the betatron coil. 

Experiment 2: FRC formation utilizing the combined beam/betatron formation 
technique. 

FRC formation was successfully demonstrated utilizing the combined beanL^betatron 
formation technique. The combined beam/betatron formation technique was performed 
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experimentally in a chamber 1 m in diameter and 1.5 m in length using an externally applied 
magnetic field of up to 500 G, a magnetic field from the betatron flux coil 320 of up to 5 kG, 
and a vacuum of 1 .2 x 10" 5 torr. In the experiment, the background plasma had a density of 
10 13 cm" 3 and the ion beam was a neutralized Hydrogen beam having a density of 1 .2 x 10 13 
cm' 3 , a velocity of 2 x 10 7 cm/s, and a pulse length of around 20 us (at half height). Field 
reversal was observed. 

Experiment 3: FRC formation utilizing the betatron formation technique. 

FRC formation utilizing the betatron formation technique was successfully 
demonstrated at the following parameter levels: 

• Vacuum chamber dimensions: about 1 m diameter, 1.5 m length. 

• Betatron coil radius of 10 cm. 

• Plasma orbit radius of 20 cm. 

• Mean external magnetic field produced in the vacuum chamber was up to 100 Gauss, 
with a ramp-up period of 1 50 ^is and a mirror ratio of 2 to 1 . (Source: Outer coils and 
betatron coils). 

• The background plasma (substantially Hydrogen gas) was characterized by a mean 
density of about 10 13 cm" 3 , kinetic temperature of less than 10 eV. 

• The lifetime of the configuration was limited by the total energy stored in the 
experiment and generally was around 30 jis. 


DOCSOCl: 147737. 1 
703538^039 K2R 


70 


PATENT 
703538.4039 

The experiments proceeded by first injecting a background plasma layer by two sets 
of coaxial cable guns mounted in a circular fashion inside the chamber. Each collection of 8 
guns was mounted on one of the two mirror coil assemblies. The guns were azimuthally 
spaced in an equidistant fashion and offset relative to the other set. This arrangement 
allowed for the guns to be fired simultaneously and thereby created an annular plasma layer. 

Upon establishment of this layer, the betatron flux coil was energized. Rising current 
in the betatron coil windings caused an increase in flux inside the coil, which gave rise to an 
azimuthal electric field curling around the betatron coil. Quick ramp-up and high current in 
the betatron flux coil produced a strong electric field, which accelerated the annular plasma 
layer and thereby induced a sizeable current. Sufficiently strong plasma current produced a 
magnetic self-field that altered the externally supplied field and caused the creation of the 
field reversed configuration. Detailed measurements with B-dot loops identified the extent, 
strength and duration of the FRC. 

An example of typical data is shown by the traces of B-dot probe signals in Fig. 30. 
The data curve A represents the absolute strength of the axial component of the magnetic 
field at the axial mid-plane (75 cm from either end plate) of the experimental chamber and at 
a radial position of 15 cm. The data curve B represents the absolute strength of the axial 
component of the magnetic field at the chamber axial mid-plane and at a radial position of 30 
cm. The curve A data set, therefore, indicates magnetic field strength inside of the fuel 
plasma layer (between betatron coil and plasma) while the curve B data set depicts the 
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magnetic field strength outside of the fuel plasma layer. The data clearly indicates that the 
inner magnetic field reverses orientation (is negative) between about 23 and 47 ^s, while the 
outer field stays positive, i.e., does not reverse orientation. The time of reversal is limited by 
the ramp-up of current in the betatron coil. Once peak current is reached in the betatron coil, 
the induced current in the fuel plasma layer starts to decrease and the FRC rapidly decays. 
Up to now the lifetime of the FRC is limited by the energy that can be stored in the 
experiment. As with the injection and trapping experiments, the system can be upgraded to 
provide longer FRG lifetime and acceleration to reactor relevant parameters. 

Overall, this technique not only produces a compact FRC, but it is also robust and 
straightforward to implement. Most importantly, the base FRC created by this method can be 
easily accelerated to any desired level of rotational energy and magnetic field strength. This 
is crucial for fusion applications and classical confinement of high-energy fuel beams. 

Experiment 4: FRC formation utilizing the betatron formation technique. 

An attempt to form an FRC utilizing the betatron formation technique has been 
performed experimentally in a chamber 1 m in diameter and 1.5 m in length using an 
externally applied magnetic field of up to 500 G, a magnetic field from the betatron flux coil 
320 of up to 5 kG, and a vacuum of 5 x lCf 6 torr. In the experiment, the background plasma 
comprised substantially Hydrogen with of a density of 10 13 cm" 3 and a lifetime of about 40 ^is. 
Field reversal was observed. 
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Fusion 

Significantly, these two techniques for forming a FRC inside of a containment system 
300 described above, or the like, can result in plasmas having properties suitable for causing 
nuclear fusion therein. More particularly, the FRC formed by these methods can be 
5 accelerated to any desired level of rotational energy and magnetic field strength. This is 

crucial for fusion applications and classical confinement of high-energy fuel beams. In the 
confinement system 300, therefore, it becomes possible to trap and confine high-energy 
plasma beams for sufficient periods of time to cause a fusion reaction therewith. 

To accommodate fusion, the FRC formed by these methods is preferably accelerated 
10 to appropriate levels of rotational energy and magnetic field strength by betatron acceleration. 
Fusion, however, tends to require a particular set of physical conditions for any reaction to 
take place. In addition, to achieve efficient bum-up of the fuel and obtain a positive energy 

r 

balance, the fuel has to be kept in this state substantially unchanged for prolonged periods of 
time. This is important, as high kinetic temperature and/or energy characterize a fusion 

15 relevant state. Creation of this state, therefore, requires sizeable input of energy, which can 

only be recovered if most of the fuel undergoes fusion. As a consequence, the confinement 
time of the fuel has to be longer than its burn time. This leads to a positive energy balance 
and consequently net energy output. 

A significant advantage of the present invention is that the confinement system and 

20 plasma described herein are capable of long confinement times, i.e., confinement times that 
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exceed fuel burn times. A typical state for fusion is, thus, characterized by the following 
physical conditions (which tend to vary based on fuel and operating mode): 

Average ion temperature: in a range of about 30 to 230 keV and preferably in a range 
of about 80 keV to 230 keV 
5 Average electron temperature: in a range of about 30 to 100 keV and preferably in a 

range of about 80 to 100 keV 

Coherent energy of the fuel beams (injected ion beams and circulating plasma beam): 
in a range of about 100 keV to 3.3 MeV and preferably in a range of about 300 keV to 3.3 
MeV. 

10 Total magnetic field: in a range of about 47.5 to 120 kG and preferably in a range of 

about 95 to 120 kG (with the externally applied field in a range of about 2.5 to 15 kG and 
preferably in a range of about 5 to 15 kG). 

Classical Confinement time: greater than the fuel burn time and preferably in a range 
of about 10 to 100 seconds. 
1 5 Fuel ion density: in a range of about 1 0 14 to less than 1 0 16 cm" 3 and preferably in a 

range of about 10 14 to 10 15 cm" 3 . 

Total Fusion Power: preferably in a range of about 50 to 450 kW/cm (power per cm 
of chamber length) 

To accommodate the fusion state illustrated above, the FRC is preferably accelerated 
20 to a level of coherent rotational energy preferably in a range of about 100 keV to 3.3 MeV, 
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and more preferably in a range of about 300 keV to 3.3 MeV, and a level of magnetic field 
strength preferably in a range of about 45 to 120 kG, and more preferably in a range of about 
90 to 1 1 5 kG. At these levels, high energy ion beams can be injected into the FRC and 
trapped to form a plasma beam layer wherein the plasma beam ions are magnetically 
confined and the plasma beam electrons are electrostatically confined. 

Preferably, the electron temperature is kept as low as practically possible to reduce the 
amount of bremsstrahlung radiation, which can, otherwise, lead to radiative energy losses. 
The electrostatic energy well of the present invention provides an effective means of 
accomplishing this. 

The ion temperature is preferably kept at a level that provides for efficient burn-up 
since the fusion cross-section is a function of ion temperature. High direct energy of the fuel 
ion beams is essential to provide classical transport as discussed in this application. It also 
minimizes the effects of instabilities on the fuel plasma. The magnetic field is consistent with 
the beam rotation energy. It is partially created by the plasma beam (self-field) and in turn 
provides the support and force to keep the plasma beam on the desired orbit. 
Fusion Products 

The fusion products are born predominantly near the null surface from where they 
emerge by diffusion towards the separatrix 84 (see Fig. 8). This is due to collisions with 
electrons (as collisions with ions do not change the center of mass and therefore do not cause 
them to change field lines). Because of their high kinetic energy (product ions have much 
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higher energy than the fuel ions), the fusion products can readily cross the separatrix 84. 
Once they are beyond the separatrix 84, they can leave along the open field lines 80 provided 
that they experience scattering from ion-ion collisions. Although this collisional process 
does not lead to diffusion, it can change the direction of the ion velocity vector such that it 
5 points parallel to the magnetic field. These open field lines 80 connect the FRC topology of 

the core with the uniform applied field provided outside the FRC topology. Product ions 
emerge on different field lines, which they follow with a distribution of energies; 
advantageously in the form of a rotating annular beam. In the strong magnetic fields found 
outside the separatrix 84 (typically around 100 kG), the product ions have an associated 

10 distribution of gyro-radii that varies from a minimum value of about 1 cm to a maximum of 

around 3 cm for the most energetic product ions. 

Initially the product ions have longitudinal as well as rotational energy characterized 
by l / 2 M(v par ) 2 and K M^) 2 . Vperp is the azimuthal velocity associated with rotation around 
a field line as the orbital center. Since the field lines spread out after leaving the vicinity of 

15 the FRC topology, the rotational energy tends to decrease while the total energy remains 

constant. This is a consequence of the adiabatic invariance of the magnetic moment of the 
product ions. It is well known in the art that charged particles orbiting in a magnetic field 
have a magnetic moment associated with their motion. In the case of particles moving along 
a slow changing magnetic field, there also exists an adiabatic invariant of the motion 

20 described by l A M(Vperp) 2 /B. The product ions orbiting around their respective field lines 
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have a magnetic moment and such an adiabatic invariant associated with their motion. Since 
B decreases by a factor of about 10 (indicated by the spreading of the field lines), it follows 
that Vperp will likewise decrease by about 3.2. Thus, by the time the product ions arrive at the 
uniform field region their rotational energy would be less than 5% of their total energy; in 
other words almost all the energy is in the longitudinal component. 

While the invention is susceptible to various modifications and alternative forms, a 
specific example thereof has been shown in the drawings and is herein described in detail. It 
should be understood, however, that the invention is not to be limited to the particular form 
disclosed, but to the contrary, the invention is to cover all modifications, equivalents, and 
alternatives falling within the spirit and scope of the appended claims. 
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